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We study numerically the photon emission from a semiconductor microcavity containing N > 2 
quantum wells under the influence of a periodic external forcing. The emission is determined by the 
interplay between external forcing and internal interaction between the wells. While the external 
forcing synchronizes the periodic motion, the internal interaction destroys it. The nonlinear term 
of the Hamiltonian supports the synchronization. The numerical results show a jump of the photon 
intensity to very large values at a certain critical value of the external forcing when the number of 
quantum wells is not too large. We discuss the dynamics of the system across this transition. 
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I. INTRODUCTION 

Over the past five decades nonlinear equations have 
been used extensively for a detailed investigation of the 
optical properties of semiconductors [ll-Q because of 
their potential application in opto-electronic devices [3~ 
Q. In semiconductor nanostructures like quantum wells 
and quantum dots [lol [ll| the coupling between light and 
matter is enhanced. By this, it may produce more pro- 
nounced nonlinear and quantum effects such as a modifi- 
cation of the quantum statistical properties of the emit- 
ted light as well as bistability and multistability. These 
effects were theoretically predicted and experimentally 
observed by several groups (l^ - l2]| . 

In natural systems as well as in experimental realiza- 
tions of artificial systems, the presence of an external 
forcing is unavoidable. An external forcing could be ei- 
ther noise (from the surrounding or inherent within the 
experiment) or caused by a deterministic perturbation. 
Sometimes, the external forcing is useful for practical ap- 
plications but sometimes it entails some degradation of 
the desired system behavior. Examples of important and 
useful results obtained by means of external perturbation 
are stochastic resonances [13: |2^, chaos control [25|. 
strange nonchaotic dynamics [26l - l28j etc. One aspect to 
explore the nonlinear behavior is to scan the parameters 
space and to observe how the dynamical complexity de- 
pends on the parameters j29l - l32| . 

The interaction between nonlinear systems gives rise 
to new phenomena such as synchronization, hysteresis, 
phase locking, phase shifting, phase-flip, riddling, ampli- 
tude death etc. [33l - l40| . Recently the coupled nonlinear 
dynamical systems have been extensively studied from 
both the theoretical and experimental points of view in a 
variety of contexts in the physical, biological, and social 
sciences etc. [H, 113] . 

In our previous work (4]| we explored the dynamics of 



the field intensities in the high excitation regime inside a 
semiconductor microcavity containing one quantum well. 
We observed periodic-doubling, quasiperiodic and direct 
route to chaos as the forcing strength is changed. These 
results show various types of dynamics depending on the 
forcing strength. Furthermore we observed also coexist- 
ing periodic and chaotic attractors with riddled basin. 

In the present paper we consider a network of quan- 
tum wells as shown schematically in Fig. 1 : quantum 
wells are inside a semiconductor microcavity (schemati- 
cally represented by two Bragg mirrors M) in the pres- 
ence of an external forcing. The aim is to analyze the 
dynamical behavior of the intra-cavity photonic intensity 
and the intensity of the fluorescent light in the presence 
of periodic signals. We observed a jump in the intensities 
of photon and exciton emission at appropriate values of 
the forcing strength. The maximum of the photon inten- 
sity is obtained for an optimal number of quantum wells 
which should not be too large. Across this transition we 
see either periodic or chaotic motion depending upon the 
exciton resonance frequency of the individual quantum 
wells. 

The paper is organized as follows. In Sec. II, we re- 
view the derivation of exciton-phonon interaction in the 
presence on A^ quantum wells. This is followed by Sec. 
Ill where the results for several quantum wells are con- 
sidered. The results are discussed and some conclusions 
are drawn in Sec. IV. 



II. MODEL 

We consider N quantum wells inside a semiconductor 
microcavity, see Fig. 1. Each quantum well is localized 
in a position, which corresponds to the maximum of the 
electromagnetic field inside the microcavity. This system 
is similar to the one described in detail in [42] . Neglecting 
the effects of the spins, the interaction with the phonons 
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FIG. 1: The schematic model for microcavity with A''- 
quantum wells. Symbols are discussed in text. 



and the excitonic saturation, we can write the Hamilto- 
nian describing the system with N quantum wells inside 
the cavity and pumped with a laser with amplitude e and 
frequency wl as |43l - l49| : 
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(1) 



The first two terms of the Hamiltonian correspond to 
the proper energies of photon and quantum well-excitons, 
where a, bj are respectively the annihilation operators of 
a photonic and excitonic modes verifying: [a, a"*"] = 1 and 
\bj, bj''\ = 1. ujph is the photonic mode frequency and ujj 
is the excitonic mode frequency for the j*'* quantum well. 
The third term corresponds to the pump energy. The 
first part of the fourth term represents the exciton-photon 
coupling with a coupling constant g^. The last part of 
the fourth term describes the excitonic nonlincarity for 
each quantum well with the coefficients a^. 

By considering that the fluctuations are weak com- 
pared to the average values, the evolution of the mean 
field operators in the interaction picture can be written 
as: 



N 



dr 

d{bj 



e(0 + XI 9] {b]) - K (a) - iAa (a) 

3 = 1 



(2a) 



-9j («) - I (b,) (6,) - 2^a, (5+) (6,) (b,) 

(2b) 

where j ~ I, N , r is a dimcnsionlcss time normalized 

to Tr, = tS- 
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7j, K are the dimensionless decay rates of the excitons 
and the cavity photon: 



(4) 



the nonlinear coupling constant aj and the coupling gj 
are normalized to — : 



9j = 9]Tc, 



(5) 



and Aq, Aj are the dimensionless detuning 

Aj (Wj-Wl)Tc. 



(6) 



According to these equations, the coupling among exci- 
tonic modes arises from their common interaction with 
the photonic mode. It represents a global mean field cou- 
pling of excitonic modes and is given by the sum term 
(over j) in Moreover, the evolution of the excitonic 
modes contains a nonlinear term (last term in ([2]). It is 
worth to mention that the nonlinear parameters aj can 
be scaled to 1 by redefining < bj""^ >=< bj > /^A^j and 
guew _ ^ Voij- This system has 2(A^ + l)-dimensions 
in the presence of forcing . For the numerical simula- 
tion we use RK4 integrator [s^l- We consider the step 
size At = 27r/5000 for integration. The dynamical stud- 
ies are studied after removing initial 10'' data points as 
transients. We explore here the photon and exciton inten- 
sities la = (a+) (a) and = (5^) (bj) inside the cavity. 
As the fluorescent light is proportional to the mean num- 
ber of excitons (Ij ) , we are also exploring the fluorescent 
light dynamics. 



III. NUMERICAL RESULTS 

The system properties are determined by different 
physical parameters. In our calculations, we fix some 
of them in order to have the possibility to determine the 
influence of the remaining ones. In all our calculations, 
the normalized parameters are fixed to gj = g = 1.5 
and K = 0.12, 7 = 7j_= 0.015 which correspond to the 
experimental values [I3 in units of the inverse of the 
round-trip in the microcavity. Furthermore, we concen- 
trate our analysis to the case of a strong pump field where 
the non-linear phenomena are expected to influence the 
dynamics. We choose the normalized amplitude of the 
laser pump as e = 200 and the normalized excitonic in- 



teraction coefficient as a,- 



0.00001. At this set of 



parameters the equations describing the dynamics of the 
system have stable fixed point solutions. 

The parameter e characterizes the influence of external 
forcing on the dynamics of the system. We use it, in the 
present study, in order to receive information on the effect 
of the external forcing. Wc restrict our analysis to a de- 
terministic periodic forcing in e, i.e. e — >■ e(l -I- / cos{Qt)) 
where Q is the perturbation frequency while / is the 
strength of the forcing. The last value is considered as 
bifurcation parameter. As to the first value, we consider 
only the case i7 = 1 in the present work. This value of 
51 corresponds to 1.5 Thz physical frequency. The Thz- 
sources arc recently realized IHJ ■ 

As we change the magnitude of /, various dynamical 
motions are possible. We have discussed this result re- 
cently |4l| for the case of one quantum well and different 
values of the detuning Aq and Ai. In the present pa- 
per we consider the effect of the forcing strength /, of 
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the detuning Aj and of the number A'' of quantum wells 
at fixed — (the cavity is resonant with the pump 
laser) . 



A. Two quantum wells 

In order to show the results in details, we first consider 
two quantum wells. We fix Ai = —g and vary A2 G 
[-9,9]- 

Shown in Figs. [Ua,b) are the photon intensities la 
in the parameter space A2 — /, respectively, without 
{aj = 0) and with {aj = 0.00001) a nonlinear term in the 
Hamiltonian ([1]) . Fig. HJa) shows that la varies smoothly 
as a function of the forcing / except near to the value of 
A2 = —0.5 when the nonlinear term is not considered. 
The results obtained with inclusion of the nonlinear term 
show another behavior. Fig. ^h). Here, the intensity in- 
creases smoothly as a function of / for all but a certain 
critical value of /. At this critical value, la jumps to 
much higher values. 

The details of the behavior of la are shown in Fig. [3] 
where the left panel is drawn for fixed Ai = A2 = g while 
the right panel corresponds to Ai = —g and A2 = g. 
These two cases correspond to identical and mismatched 
quantum wells, respectively. The mismatch in the latter 
case is equal to the Rabi frequency of the single quantum 
well which is |Ai — A2I = 2g. 

The results shown in Figs. [2Ja,b) are obtained with 
aj = 0, i.e. with a vanishing nonlinear term in the Hamil- 
tonian H, Eq. ([1]). They show an overall smooth increase 
of the intensity with increasing forcing strength /. The 
results in Figs. [2Jc,d) are obtained with a nonvanishing 
nonlinear term in ([T]), i.e. with aj = a = 0.00001. In 
this case, la does not increase everywhere smoothly with 
increasing forcing /, see also Fig. [2][b). The nonlinearity 
causes a substantial jump in the intensity around / ^ 2. 
This jump occurs independently of the detuning A2. In 
Figs. [3Ic,d) the jump is marked by an arrow. 

In order to see the dynamical behavior of the system 
across this transition we plot a few of the largest Lya- 
punov exponents (LEs) in Figs. [3][e,f) for the nonlinear 
cases considered in Figs, ^c) and (d), respectively. The 
LEs arc calculated according to [53 • In the case of iden- 
tical quantum wells (left panel) the dynamics is always 
periodic and all but one Lyapunov exponent (which dis- 
tinguish the type of dynamics, see Ref. [1^) are neg- 
ative. One of the LEs is zero (dotted line). The third 
LE (dashed line) is negative but jumps at the same value 
/ = 1.3 at which the intensity jumps (shown by an arrow 
in Fig. |3l[c)). The trajectories in the phase space across 
this transition arc shown for two different values of the 
forcing strength / in Fig. |3|g) . The inner black solid line 
corresponds to / = 1 while the outer red-dashed line is 
for / = 2. These results are confirmed by the correspond- 
ing Poincare sections [2^ which are taken at Re{a) = 0, 
see Fig. ^i). There are two single points corresponding 
to / = 1 and 2. 




FIG. 2: The photon intensity la in the parameter space A2 — / 
at the fixed value Ai = —g for two quantum wells without 
(a) and with (b) nonlinear term in the Hamiltonian Eq. ([T]). 
The contour line in (b) is drawn at la ~ 2 x 10^. 



It should be noted here that both the photon intensity 
and the exciton intensity jump at the same values of /. 
The photon intensity is, however, much larger than the 
exciton intensity, see Fig. El^c). The variation as a func- 
tion of time is shown in the insets of Fig. EJe) which are 
calculated with the parameters of Fig. [SJg). 

In the mismatched case (right panel) we see a behavior 
of the intensities / as a function of the forcing strength 
/ which is similar to that discussed above for the case 
with identical wells. Due to the mismatching, the crit- 
ical values differ from those of the left panel. However 
the exciton intensities are much smaller than the pho- 
ton intensities also in this case. The dynamics across the 
transition is shown in Fig. ^i). The spectrum of the 
LEs shows the following characteristic features : below 
the transition, the motion is chaotic while it is periodic 
beyond the transition. Both, the chaotic and the periodic 
motion arc shown in Fig. El^h) with black-solid and red- 
dashed lines, respectively, at the forcing strengths / = 1 
and 2. These results arc confirmed by the correspond- 
ing Poincare sections taken at Re{a) = and shown in 
Fig. [31Jj). The motion is chaotic (scattered points) and 
periodic (single point), respectively. The variation of the 
intensities as a function of time is shown in the insets of 
Fig. at the parameter values / = 1 and 2. 

Comparing the results obtained for the case with two 
identical quantum wells (left panel of Fig. [3]) to those 
obtained with two mismatched wells (right panel of Fig. 
[3]), we state the following. In the first case, the intensities 
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jump at a certain critical value of the forcing strength / 
while the transition starts much below this critical value 
of / in the second case. Note that the y-axis is taken 
on logarithmic scale. Further, there is a synchronized 
periodic motion (where the Ij have identical values [s^) 
across the jump in the first case with identical wells while 
the motion changes from an unsynchronized chaotic to a 
synchronized periodic one in the second case with mis- 
matched wells. In any case, the intensities / at large / 
are much larger when the excitonic nonlinearity in the 
Hamiltonian Eq. ([T]) is taken into account, than with- 
out this term. The nonlinearity causes, obviously, the 
jump-like transition to the higher intensities at the criti- 
cal value of the forcing strength /. 



B. Chain of N quantum wells 

In order to see the effect of the nonlinearity on the pho- 
ton and exciton intensities in the case of a large number 
N of quantum wells we consider in the following a chain 
of identical quantum wells with fixed Aj = Vj. In Fig. 
Slthe calculated photon intensity la is drawn in the pa- 
rameter space : number N of quantum wells and forcing 
strength /. 

According to Fig. HJa) the intensity la decreases 
smoothly with increasing N (> 2) when the nonlinear 
term in ([T]) in not taken into account (corresponding to 
Uj ~ 0). However there is a jump in intensity from N = 1 
to N = 2. Fig. \Mjo) shows the results obtained with in- 
clusion of the nonlinear term in ([T]), i.e. with aj ^ 0. 
In this case, the intensity increases first with increasing 
iV and then jumps to lower values. This behavior re- 
peats several times. Both figures indicate further that, 
at very low forcing strength f 0, la does not change 
for any value of N. As a result, a jump in the intensity 
la appears only when the forcing strength / as well as 
the nonlinearity a do not vanish. 

When the nonlinearity in ([1} is taken into account in 
the calculations, a substantial drop in the intensity la 
appears at some values of the forcing strength / when we 
increase the number N of quantum wells. This behavior 
of la is determined obviously by the nonlinear term in ([l} . 
The details of the variation of the photon intensity la as 
well as of the exciton intensity Ij (all Ij are synchronized) 
are shown in Fig. [SJ 

Figs. [5ja, b) show numerical results obtained without 
the nonlinearity in ((U as a function of N and /, respec- 
tively. Here, the intensities vary smoothly : the photon 
intensity la as well as the exciton intensity Ij decrease 
with increasing N but increase with increasing /. The 
situation is, however, completely different when the non- 
linearity is taken into account. In this case, substantial 
jumps appear in the intensities, see Figs. [5Kc,d) and the 
further results shown in Figs. [5je-h). 

Let us first consider the variation of la as a function 
of N for fixed forcing strength / = 15. According to Fig. 
[S](c), la increases first with the number N of quantum 
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FIG. 3: (Color online) Comparison of some results for two 
identical quantum wells with those for two mismatched wells. 
The left and right panels correspond to Ai = A2 = g and 
Ai = — (7, A2 = g, respectively, (a-d): Intensities la (o), Ii 
(o) , and I2 {*) as a function of the forcing strength / without 
(a,b) and with (c,d) nonlinear term in ([T]). The synchronized 
intensities Ii (o) and I2 (*) in (a) and (c) are overlapping. 
(e,f): A few largest Lyapunov exponents as a function of the 
forcing strength /. The dotted lines represent the zero Lya- 
punov exponent. Insets in (e,f): intensities la (black-solid 
line), Ii (red-dashed line) and I2 (blue-long-dashed line) as a 
function of time corresponding to (g,h). (g,h): Trajectories in 
the phase space Im{a) — Re{a) at / = 1 (inner solid line) and 
/ = 2 (outer red-dashed line) below and above, respectively, 
the transition. (i,j): The Poincare section of (g,h). The other 
parameters are the same as those in Fig. 2. 



wells. However, it reaches a maximum value at iV = 6 
where it decreases suddenly and drops to very low values 
(even lower than in the single quantum well). The details 
of the dynamics across this transition are shown in Fig. 
[SJc). Here, trajectories at iV = 4 (outer red-dashed line) 
and 8 (inner black-solid line) are drawn in the phase space 
Re{a) — Im{a). These trajectories show periodic motions 
what is confirmed by the Poincare section shown in Fig. 
\51e) ■^ith points corresponding to TV = 4 and 8. 

Fig. [SKd) shows the intensity la as a function of the 
forcing strength / for a fixed number TV = 6 of quantum 
wells. Similar to the case with two identical wells (Fig. 
^c)) the intensity la increases smoothly up to a certain 
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FIG. 4: (Color online) The photon intensity in the pa- 
rameter space f — N for A'' identical quantum wells at fixed 
Aj = g, Vj without (a) and with (b) nonlinear term in ([T]). 



value of / (marked by an arrow in the figure) where it 
jumps to a mucli larger value. The corresponding dy- 
namics in phase space Re{a) — Im{a) is shown in Fig. 
Eljf) for / = 10 (inner black solid line) and / = 20 (outer 
red-dashed line). Both motions are periodic. This result 
is confirmed by the Poincare section, see Fig. [Hh) where 
the two corresponding points / = ID and 20 arc shown. 

The exciton intensities Ij show a similar behavior as 
the photon intensity la in all cases. They are however 
smaller than /q. 

The results shown in Fig. [5] indicate that the dynam- 
ics of the system consisting of N forced coupled semicon- 
ductor microwave cavities is determined by two opposite 
tendencies. On the one hand, the interaction between the 
quantum wells prevents a synchronized periodic motion. 
On the other hand, however, the nonlinear terms in the 
Hamiltonian ([T]) support the synchronized motion. As a 
result of these two tendencies, it is possible to increase 
the photon intensity la substantially when the number 
of quantum wells is not too large. 
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FIG. 5: (Color online) Some results for A'^ identical quantum 
wells as a function of the number A'^ (at fixed /, left panel) 
and of the forcing strength / (at fixed A'', right panel), (a-d): 
Intensities la (o) and Ij (★) as a function of A' at fixed forcing 
strength / = 15 (a,c) and as a function of / at fixed number 
Af = 6 of quantum wells (b,d). The results are obtained, 
respectively, without (a,b) and with (c,d) nonlinear term in 
(IT|. (e,f): Trajectories in phase space Im{a) — Re{a) at (e) 
N — 4 (outer red-dashed line) and A'' = 8 (inner black solid 
line) for fixed / = 15; and at (f) / — 10 (inner solid line) and 
/ = 20 (outer red-dashed line) for fixed A'' = 6. (g,h) The 
Poincare section of (e,f). The other parameters are the same 
as those in Fig. 4. 



IV. DISCUSSION AND CONCLUSIONS 

In Sect. Illli we showed the calculated intensities of 
photon and exciton emission from a semiconductor mi- 
crocavity containing of N quantum wells under the in- 
fluence of an external forcing /. The intensity of the 
photon and exciton emission, la and Ij respectively, is 
determined by the interplay of external forcing and in- 
ternal interaction between the single quantum wells of 
the microcavity. External forcing synchronizes the peri- 
odic motions and causes, by this means, an enhancement 
of the intensities. The internal interaction, however, dis- 
turbs the synchronized motion and leads to a reduction 
of the intensities of photon as well of exciton emission. 
The interplay between external and internal interaction 



is described well by the Hamiltonian 

In the case of A^ = 2, the destructive role of the inter- 
nal interaction is small. The intensities la and Ij increase 
with / starting at a certain small value of /. This holds 
true for the case of two identical wells as well as for two 
mismatched wells. In both cases, the intensities jump 
to much higher values at a certain critical value of the 
forcing strength /. In the first case, the dynamics of the 
system is a synchronized periodic motion below as well as 
beyond the jump. In the second case, however, the mo- 
tion is unsynchronized chaotic below the critical value of 
/ where the jump occurs, and changes to a synchronized 
periodic motion for / values beyond the jump. 

In the case with A^ wells, the destructive role of the in- 
ternal interaction between the wells can directly be seen. 
Without the nonlinear term in the Hamiltonian ([T]), the 
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intensities decrease with increasing TV (and fixed /) but 
increase with / (and fixed N), see Figs. [Ha,b). With 
nonlinear term in ([T}, however, the intensities increase 
first smoothly up to a certain critical value of N (and 
fixed /), and a critical value of / (and fixed N), respec- 
tively, see Figs. [5Ic,d). Beyond these values, the inten- 
sities decrease abrupt in the first case (Fig. EKc)) while 
they jump to higher values in the second case (Fig. EKd)). 
In both cases, the intensity la jumps by several orders of 
magnitude at the critical points. 

These results illustrate very nicely the interplay of in- 
ternal interaction between the wells and external forc- 
ing. The jump in the intensities is a coherent effect that 
does not occur in a single quantum well, see the results 
obtained earlier [4ll |. The jump in the intensities with 
a fixed (not too large number N of wells) at a critical 
value of / is similar to that observed in the case with 



two wells. In any case, the dynamics across the jump de- 
pends on the internal parameters of the quantum wells, 
and the nonlinear term in the Hamiltonian ([l} plays an 
important role. 

The enhancement of the intensity of photons emitted 
from a microcavity is of great interest for applications. 
According to the results of the present paper it is possi- 
ble to manipulate microcavities in such a way that the 
intensity of emitted light is very large. 
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and acknowledge the financial support and the hospital- 
ity of the MPIPKS. 
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